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We study the behavior of the AsqTad quark propagator in Landau gauge on quenched SU(2) gauge 
configurations under the removal of center vortices. In contrast to recent results in SU(3), we clearly 
see the infrared enhancement of the mass function disappear if center vortices are removed, a sign 
of the intimate relation between center vortices and chiral symmetry breaking in SU(2) gauge-field 
theory. These results provide a benchmark with which to interpret the SU(3) results. In addition, 
we consider vortex-only configurations. On those, the quark dressing function behaves roughly as 
on the full configurations, and the mass function picks up an almost linear momentum dependence. 
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I. INTRODUCTION 

Dynamical breaking of chiral symmetry (D^SB) is an 
essential nonperturbative property of quantum chromo- 
dynamics (QCD) which cannot be accounted for within 
perturbation theory at any order. Only nonperturbative 
approaches, such, as those provided by lattice QCD simu- 
lations or studies of the Dyson-Schwinger equations, can 
be used to explore this phenomenon. 

The other characteristic nonperturbative phenomenon 
of QCD is confinement: the fact that colored states 
are never observed. It is tempting to speculate that 
these two phenomena might be driven by the same ba- 
sic mechanism, an idea supported by finite-temperature 
studies where the deconfinement and chiral restoration 
transitions are observed to occur at coincident temper- 
atures Moreover, it was found that the low- lying 
modes of the quark operator not only bear witness on 
spontaneous chiral symmetry breaking but also on con- 
finement 

One leading candidate for such a mechanism is the cen- 
ter vortex. Center vortices have been studied in lattice 
QCD simulations for more than a decade. The recovery 
of the string tension from "vortex-only" SU(2) gauge con- 
figurations [i.e., Z2 projected from SU(2)] is well known, 
as is the recovery of the chiral condensate [1, 0, H, @| 
(see also more recent Refs. 0, H, Q). In SU(3), how- 
ever, the situation is far less compelling. For example, 
in Refs. [HI numerical evidence has been given that 
in SU(3) mass generation remains intact after removing 
center vortices, whereas the string tension vanishes as 
expected. This came as a surprise and immediately sug- 
gested a corresponding "bench-mark" study for the case 
of SU(2) which is reported here. 

We use the quark propagator as a probe of D^SB. The 
Dirac scalar part of the propagator, related at large mo- 
menta to the perturbative running mass, is enhanced at 
low momenta, even in the chiral limit Il3|: a demon- 



stration of DxSB. We explicitly establish the relation 
between center vortices and D^SB by investigating the 
quark propagator under the removal of center vortices. 
Specifically, we will provide numerical evidence that dy- 
namical mass generation disappears if those vortices are 
removed, and surprisingly, much of it resides in the 
vortex-only part. Additionally, we present an improved 
method for generating the SU(2) maximal-center-gauge- 
projected configurations. 

The paper is organized as follows. In an attempt to 
make it self-contained we briefiy introduce the quark 
propagator and its realization on the lattices based on 
the AsqTad quark action. This is followed by a specifi- 
cation on how wc gauge-fixed and identified center vor- 
tices in our gauge configurations. Then, results for the 
mass and quark dressing function are compared on full, 
vortex-removed and vortex-only configurations. A sum- 
mary concludes the paper. 



II. THE QUARK PROPAGATOR 

The quark propagator is gauge dependent. In co- 
variant gauges in Euclidean momentum space it can be 
parametrized in the general form 



Sip') 



Zip') 



(1) 



where M is the running mass and Z the quark dressing 
function. S can be calculated in regularized theories, as 
for example here in a lattice regularization where the lat- 
tice spacing a makes all expressions finite. At sufficiently 
small a, i.e., if scaling violations due to finite lattice spac- 
ings arc negligible, the bare quark propagator, S, is re- 
lated to the renormalized propagator via multiplicative 
renormalization: 

SRip'; IJ., gRifJ.),mR{fi)) = Z2S{p'^;a,go{a),77io{a)). 
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To ensure multiplicative renormalizability, all the depen- 
dence of Sn on the renormalization point /i is contained 
in the renormalized quark dressing function, Z^. M does 
not depend on /i. A renormalization condition fixes Z2, 
the renormalization constant of the quark fields. Lattice 
calculations often use MOM schemes to fix renormaliza- 
tion constants. In MOM schemes, Z2 is fixed by requiring 
Sr to be of the form of a free propagator at the renor- 
malization point = This sets 



= 1 and = mnin) 



(2) 



where the latter denotes the renormalized mass at /i^. 

Calculation of the quark propagator S proceeds like 
any correlation function in a lattice Monte Carlo (MC) 
calculation once the gauge has been fixed. For the gauge 
we used the ever popular Landau gauge. It is straight- 
forward to implement on the lattice and allows for an 
easy comparison to other studies. The Landau-gauge 
quark propagator has been studied widely in SU(3) gauge 
theory using Wilson-clover, staggered-type and overlap 
actions in quenched and unquenched simulations (see, 
e.g, Refs. [H [H, M, M, m [13, E, E Hi), it has 
been shown that the quark propagator obtained with 
the AsqTad-improved staggered quark action possesses 
good symmetry properties and is well behaved at large 
momenta (see, e.g., Ref. [U). The AsqTad action was 
therefore a natural choice for this study. 



III. DETAILS OF THE CALCULATION 

Configurations were generated on a 16'^ x 32 lattice 
using a tadpole-improved Wilson gauge action with an 
inverse coupling constant /3 = 1.35. Around 120 configu- 
rations were used. 



A. Identifying center vortices 

The observation that the long-range static potential 
only depends on the center charge (also called 'W-ality" ) 
of the quark representation led to the expectation that 
the center subgroup Z^ of SU(A^) plays a crucial role 
for quark confinement. Center vortices emerge from the 
corresponding Zn gauge theory: they form closed world- 
sheets in four space-time dimensions carrying flux which 
takes values in the center of the SU(A^) group. Early at- 
tempts to define this subgroup by the projection S\J{N) 
Ztq failed in the sense that the arising vortex mat- 
ter did not have meaningful properties in the continuum 
limit of vanishing lattice spacing a. In the pioneering 
works [13, a two step process was proposed to define 
the links Zfj_(x) spanning the Zn gauge theory: 



(i) ^|Tr[/^«(x)|^ ^ max 



X,fJ. 



(ii) D\cTr(ui',{x)Zl{. 
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FIG. 1: Scaling of the planar vortex area density p in units 
of the string tension a as function of the lattice spacing a. 



Step (i) is the difficult part involving gauge fixing and 
the interference of Gribov ambiguities (see comments be- 
low). The projection step (ii) operates locally and can be 
implemented straightforwardly. Using a standard itera- 
tion overrelaxation procedure [1^ for step (i) finally de- 
fines vortex matter with a sensible phenomenology in the 
continuum limit: for SU(2), it was observed that the so- 
called "vortex-only" configurations, defined by the center 
links Zi\[{x) reproduce a good deal of the string tension 
while "vortex-removed" configurations, spanned by the 
links 



(5) 



do not support a linear rising static potential at large 
distances [22|, . It was put forward in [2^ that this 

heuristically defined vortex texture has sensible proper- 
ties in the continuum limit. It was subsequently dis- 
covered that these vortices provide an intriguing picture 
of the deconfinement phase transition at finite temper- 
atures [2^. They even admit detailed insights in the 
critical phenomenon [27l . [2^ in complete agreement with 
studies of the free energy of gauge-invariant center vor- 
tices over this transition [29j . 

We point out that present day algorithms [s^l are only 
capable of obtaining one of the many local maxima of the 
gauge-fixing condition ^ and that the vortex properties 
do depend on the set of maxima [sil . [s^ . [33j . Localizing 
the global maximum of Eq. ^ would remove this ambi- 
guity, but it is not clear whether the vortex matter arising 
from the global maximum would be of any phenomeno- 
logical value. In fact, an anticorrelation was found [33 |: 
the larger the gauge-fixing functional ([3]), the lesser the 
string tension obtained from vortex-only configurations. 
Although we have not yet found a concise mathematical 
description, the vortex matter best for phenomenologi- 



3 



cal studies seems to arise from an average over Gribov 
copies. 

Note also that the relation between vortices and quark 
confinement is far less striking for the gauge group SU(3): 
despite many attempts, the string tension arising from 
vortex-only configurations i s sy stematically smaller than 
the full string tension [^,[33,|3i|- The difference in quality 
between the SU(2) and SU(3) vortex picture indicates 
that we are still missing a point for SU(3) and partially 
motivated the present investigation of the SU(2) quark 
propagator. 

Because of the above ambiguity, we are going to 
present the details of the gauge-fixing procedure which 
gives rise to an intriguing vortex phenomenology. For 
SU(2), the gauge-fixing matrix g{x) £ SU(2) can be 
viewed as a 4-dimensional unit-vector in Euclidean color 
space: 

g{x) = goix) + ^fg{x) , G{x) = ■'^^(^^^ ) , 
where 

gl{x)+f{x) = l. 



Maximizing the gauge-fixing functional ([3]) is an itera- 
tive procedure: choosing a particular site Xq and setting 
g{x 7^ xq) = 1, the functional is locally maximized, and 
the links are updated accordingly: — > J7^. Subse- 
quently, all sites are visited and many sweeps through 
the lattice are performed until the gauge-fixing action 
does not change anymore within the required precision. 

For the local update (?(a;o), let sgx denote that part 
of the gauge-fixing functional ([3|) which is affected by a 
change of g{xQ): 



Sfi, = G^(xo) M G(xo) - A (g^{x^)G{x^) - l) , 



where A is a Lagrange multiplier and M is a real sym- 
metric 4x4 matrix given in terms of the link fields 
[Uf_i{x) = u°(a;) -I- iTu^j,{x)\. 



M{x) 



4 

E 



{ul{x)f + {ulix - ^l)f -ul{x)u^^{x) -f ulix ~ fi)ul{x - 

~ul{x)ul{x) + u^ix - - /i) u'^{x)ul{x) + u^ix - n)ul{x - fl) 



(6) 



We also introduce the eigenvectors and eigenvalues of the 
matrix M : 



Gfc, A*. 



with fc = 1 . . . 4 



Choosing the largest eigenvector for the gauge trans- 
formation, i.e., G{xo) = eijiax, the local increase of the 
gauge-fixing functional is maximized: 

This choice gives rise to the standard iteration procedure 
which is usually employed for MCG fixing [2^. 

Here, we depart from this standard procedure and in- 
troduce an aspect of simulated annealing: we choose 
G{xo) — Ck with a relative probability of exp{/3/Afc}, 
where f3f is an auxiliary parameter familiar from simu- 
lated annealing. For large (3f, the probability for picking 
the largest eigenvalue is high, and our method smoothly 
merges with the standard scheme. In practice, we started 
with /?/ = 0.02 and performed 25 sweeps through the lat- 
tice until we increased /3f by 0.1. The procedure stopped 
when no further increase of the gauge-fixing functional 
was achieved. 

It turns out that the vortex matter arising from this 
procedure has good phenomenological properties such as 



good scaling properties in the continuum limit. To test 
the latter aspect, we calculated the planar vortex area 
density p in units of the (measured) string tension for 
several values of the lattice spacing a using the stan- 
dard Wilson action. Figure [1] illustrates how the vortex 
density becomes independent of the lattice regulator for 
sufficiently small values of the lattice spacing. 



B. Fixing to Landau gauge 

To fix configurations to Landau gauge an overrelax- 
ation algorithm was used. This is an iterative algorithm 
that maximizes the Landau-gauge functional 



Fu[g] 



1 



EfteTr 



9xUx,t,g. 



t 

X+fl 



(7) 



by changing the gauge transformation fields g^ locally 
but keeping U fixed. The algorithm stopped if an accu- 
racy of 



maxTi- [df^A.^^d^Al J < 10 



-13 



(8) 
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was reached [4l|. Here the gauge potential is defined as 
1 



A — 



Ziago 



For the vortex-only configurations, Fourier acceleration 
was also used [36| . We applied the standard technique 
where each configuration was gauge-fixed once. 



C. The AsqTad quark propagator 

On those gauge-fixed configurations we calculate the 
quark propagator y) in coordinate space by inverting 
the AsqTad fermion matrix. After Fourier-transforming 
S{x,Q), the mass and dressing functions, M and Z, are 
extracted from S by suitable projections in Dirac space 
(see below). At this step it is important to know how the 
discrete lattice momenta 

P^i = 7rfcp/ip with k^, e (-L^/2,L^/2] 

are related to the physical momenta o?q^(p) in lattice 
units. From experience we know it is always good prac- 
tice to look at the tree-level form of the lattice propa- 
gators, and to define momenta such that the continuum 
tree-level expression is retrieved. This, known as tree- 
level correction [sj, accounts for the lowest order dis- 
cretization effects. With the AsqTad action the tree-level 
form of the quark propagator in Landau gauge reads 

4 

'S'l ^(aV) = i X! ^A^a^MbM) + 



where is as above and 

aqt^ip,,) = sin(|j^) 



(9) 



is the "kinematic momentum" (see, e.g., Ref. p^). 
The four matrices 7^ form a staggered Dirac algebra 
(Eqs. (A.6) and (A.7) of Ref. [H]). Consequently, we 
use Eq. ^ to define momenta. 

The functions M and Z are then extracted from Sl (p) 
by taking the two traces in Dirac space 



AL{q,a) 



■tr(7^g^S'L(q,a)) 



1 



BL{q,a) ^ — tr5L(q,a) 

and combining them to (here Nc = 2) 

q^a^Al{q,a) + Bl{q) 



and 



ZL{q,a) 



ML{q,a) 



Ahiq, a) 



ALiq,a) 



(10) 
(11) 

(12) 
(13) 



In the region of asymptotic scaling. Ml {q, a) becomes 
independent of a and equals the running mass. Z needs 
to be renormalized as explained above. As we are 
only interested in qualitative changes of the momentum- 
dependence of M and Z under the removal of center 
vortices, we prefer to work directly with unrenormalized 
quantities. Therefore, the presentation of the data is 
simplified by considering only bare lattice functions at 
a fixed lattice spacing. Consequently, lattice artefacts 
could not be directly assessed; however, for this value of 
the coupling, discretization errors are not expected to be 
significant. 

A cylinder cut [sj is applied to all the data to reduce 
the effects of rotational symmetry violation. 



IV. RESULTS 
A. The mass function 

In Fig. [5] we compare the mass function on our sets 
of full, vortex-only and vortex-removed configurations. 
Data was obtained for a range of bare quark masses from 
ma = 0.020 to ma = 0.100, and the two extremes are 
shown here as functions of momentum. 

Our results on the full SU(2) configurations (see Fig.[2l 
left column) show a large enhancement near zero momen- 
tum, while data drops rapidly to its expected asymptotic 
behavior at large momentum. Also, as expected, the in- 
frared enhancement is stronger for the smaller bare quark 
mass. That is, our data for M(g^) on the full configu- 
rations clearly reproduce the well-known characteristics 
expected for the mass function in both the nonperturba- 
tivc and perturbative regime. 

On the vortex- removed configurations (see Fig. [21 right 
column), the mass function is more or less flat taking 
values slightly above ma. The momentum dependence 
of M is almost linear with a bigger slope for ma = 0.1. 
Thus, the dynamical contribution to the mass function, 
which we clearly see on the full configurations, disappears 
when center vortices are removed. 

Interestingly, the mass function on the vortex-only con- 
figurations (see Fig. [21 middle column) depends quite 
strongly on momentum. Even though the signal is quite 
noisy, M(q^) grows almost linearly upon decreasing a^q^. 
That is, much of the infrared enhancement of M is con- 
tained in the vortex-only part which clearly underlines 
the importance of center vortices as IR degrees of free- 
dom. 



B. The quark dressing function 

The same comparison for the bare quark dressing func- 
tion, Z, is shown in Fig. [31 Again, we show data at 
ma ~ 0.02 and ma = 0.10 on full SU(2) configuration 
in the left column, while the middle and right column 
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FIG. 2: The mass function aM{q) versus kinematic momentum for two bare quark masses ma = 0.02 (top) and ma = 0.10 
(bottom). We show data on fuU, vortex-only and vortex-removed configurations from left to right. Data has been cylinder cut. 
The nonperturbative enhancement of the mass function at low momenta is associated with the presence of center vortices. 



displays data on vortex-only and vortex-removed config- 
urations, respectively. As expected, on the full configu- 
rations, the quark dressing function takes values around 
one at large momenta and becomes suppressed towards 
lower momenta. The smaller the quark mass the more 
pronounced the dip at lower momenta. If center vortices 
are removed, the infrared suppression disappears and Z 
is a flat function of momentum (see the right column). 
There, Z roughly stays at its tree-level value. Surpris- 
ingly, on the vortex-only configurations, Z has a similar 
momentum dependence to the full configurations. Again, 
the results are noisier, but the infrared suppression is un- 
ambiguous. 



v. CONCLUSIONS 

We have studied the Landau-gauge quark propagator 
in quenched SU(2) gauge theory under the removal of 
center vortices. Our implementation of this propagator 
is based on the AsqTad-improvcd staggered quark ac- 
tion modified to SU(2). The full propagator is found to 
strongly resemble that of the SU(3) theory. 



Our results for the mass and quark dressing functions 
unambiguously show the disappearance of D^SB when 
center vortices are removed. This is in contrast to the 
situation in SU(3) [l3, El- There, even after center- 
vortex removal dynamical mass generation survives while 
the string tension is flat. 

Additionally, we have studied the quark propagator 
on vortex-only configurations. Even though the signal 
is quite noisy, both parts of the propagator reveal a 
form qualitatively similar to the full, untouched configu- 
rations. 

In summary, our SU(2) results clearly represent a 
strong relationship between the vortex picture and spon- 
taneous breaking of chiral symmetry. We do not expect 
this to be merely a feature of staggered fermions, how- 
ever it might be valuable to repeat this study with an- 
other type of fermion, in particular after the new results 
presented in Ref . Q . More importantly, whether it is pos- 
sible to observe similar behavior in SU(3) [or ultimately 
SU(A^)] is still an open question. 
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FIG. 3: The bare quark dressing function, Z, as a function of momentum for two bare quark masses ma = 0.02 (top) and 
ma = 0.10 (bottom). The left column displays data on full configurations, whereas in the middle and right column the data 
on vortex-only and vortex-removed configurations are shown, respectively. Data has been cylinder cut. The nonperturbative 
suppression at low momenta is associated with the presence of center vortices. 
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